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Abstract 

The polarization tensor is calculated which originates from interaction 
of a photon with the electron-positron field. The effects of multiple scat- 
tering of electrons and positrons in a medium side by side with external 
fields is included. The corresponding general representation of polarization 
tensor is found in the form of two-dimensional integral. The both effects 
may be essential for propagation of high-energy photons in oriented single 
crystals. 
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1 Introduction 



As known, the propagation of electromagnetic wave in a medium is denned by 
its dielectric tensor £^{u). For relatively low frequency u (e.g. visible light) the 
dielectric tensor is defined by atomic phenomena. When the frequency of wave 
is much higher the atomic frequencies, the dielectric tensor has a form 

UJ 2 4:7171 & 2 

£ik(w) = 6ik£(u), £(u) = 1 £, u 2 = — , (1.1) 

where n e is the number density of electrons in a medium, m(e) is the electron 
mass (charge), ujq is the plasma frequency, in any medium ojq < 100 eV. So that 
for u ^> Uq an influence of atomic phenomena on propagation of electromagnetic 
wave in a medium becomes small. It should be noted that for radiation of high- 
energy particles the characteristic frequency is 7cj , 7 is the Lorentz factor. 

At very high energy the nonlinear effects of QED enter into game. These 
effects are due to interaction of photon with the electron-positron field. One 
of them is the polarization of the vacuum by a photon. In the presence of an 
external electromagnetic field the polarization of vacuum was considered first 
in the pioneer papers 0. In the strong field this effect can be essential for 
propagation of high-energy photons @, @. 

To evaluate the polarization tensor one has to consider the amplitude of pho- 
ton scattering which included the polarization operator. We use the quasiclassical 
operator method ||, ||, [H, ||. In this method the mentioned amplitude is de- 
scribed by diagram where the virtual electron-positron pair is first created by the 
initial photon with 4- momentum k(u, k) and polarization ei and then annihilate 
into final photon with 4-momentum k and polarization e 2 . This corresponds to 
use of the non-covariant perturbation theory where at high energies (u 3> m) 
the contribution of this diagram survives only. For this energy of photon this 
process occurs in a rather long time (or at a rather long distance) known as the 
lifetime of the virtual state 

where q c > m is the characteristic transverse momentum of the process, the 
system h — c — 1 is used. When the virtual electron (or positron) is moving in a 
medium it scatters on atoms and changes the velocity under influence of external 
electromagnetic field. The mean square of momentum transfer to the electron 
from a medium and an external field on the distance // is 

If = Qs + (ll = A7rZ 2 a 2 n a Ll f} L = L(q 2 ) = \n(q 2 a 2 ), 

q F = eF e l f , F e = E ± + vxH (1.3) 

where a = e 2 = 1/137, Z is the charge of nucleus, n a is the number density of 
atoms in the medium, a is the screening radius of atom, is the electric field 
strength transverse to the velocity of particle v ~ n = k/cu, H is the magnetic 
field strength. 
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In the case of small momentum transfer q/ = yq 2 <C m the influence of a 
medium and an external field is weak, in this case q c = m. At high energy it 
is possible that q c > m. In this case the characteristic value of the momentum 
transfer (giving the main contribution into the spectral probability) is defined 
by the value of qf. The self-consistency condition is 



9 9 27iu)Z 2 a 2 n a L(qf) m e K, 2 9 u „ 



If 4( 1. 



f 



F = E — n(nE) + n x H, (1.4) 

here k = \k\ is known parameter characterising the pair production process in a 
homogeneous external field F. With q c increase the lifetime of the virtual state 
(|1.2| ) decreases. 

We will use the following normalization condition for the amplitude under 
consideration 

M = 2luAu. (1.5) 

The amplitude M is the contraction of the tensor ej*e k (e^ and e^) are 
the polarization vectors of the initial and the final photons) and the polarization 
tensor Mj k - We select the basic vectors as 

F 

ei = — , e 2 = n x e x . (1.6) 
r | 

Since the tensor Mjk is invariant under the space inversion then in the selected 
basic vectors it has the diagonal form 

M ik = ^ [8jk(M n + M 22 ) + ( eii eifc - e 2j e 2k )(M u - M 22 )] (1.7) 

In absence of external field it is convenient to describe the process of photon 
scattering using the helicity polarization vector eA (A = ±1) connected with 
momentum transfer A (see Eq.(2.34) in ||). In presence of external field and 
for A = we choose the polarization vectors in the following way: 

e A = i (ei + i\e 2 ) , (e A e A ) = 1, (e A e* A ) = 0, e A x n = iAe A . (1.8) 

In terms of helicity amplitudes M ++ and the tensor Mjk and the corre- 
sponding dielectric tensor Ejk has a form 

Mj k = 5 jk k 2 ++ + (eije lk - e 2i e 2fc )A;+_, 

£ jk = 6 jk -— M jk , k 2 ++ = M ++ , k 2 + _ = M + _ (1.9) 

UJ 

The polarization tensor k 2 k is diagonal in the basic vectors ei and e 2 (|1.6|) . The 
corresponding mass squared are 

k\ = k\ x = k 2 lr+ + k 2 ^__, k\ = k 2 2 = k 2 _ + — k 2 l[ __ (1-10) 
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The probability of pair creation by a photon with polarization e is 



<(e) 



1 



Im 



eexYkf + (ee 2 ) 2 /c. 



2 ;„2 
2 



Im 



where £3 is the Stokes' parameter. For unpolarized photon one has 



(1.11) 



■— Im k 2 ,, . 

UJ 



1.12) 



2 Polarization tensor 



The polarization of a medium by a high-energy photon in the presence of an 
external field is described by the amplitudes of photon forward scattering. The 
general representation of photon scattering amplitudes without change of helic- 

ity M ++ = M and with helicity flip M + _ = M_ + was obtained in || (see 

Eq.(3.14)). In the mentioned paper || the helicity amplitudes were normalized 
by the condition ImM ++ = unr p (u), where a p {uj) is the total cross section of a 
production of electron-positron pair by a photon with energy u. In the case un- 
der consideration we are interested in corrections to the energy (mass) of photon 
and it is convenient to use the normalization ( 1.11 ) and (|1.12|) 



k 2 



-2am 2 / T 



code 



XX' 



ee' 



(2.1) 



where 

T + 
T + 

T + 

S3 



= T__, T+_ = T_ 

= (0|s a (G- 1 - G, 1 ) + s 2 p (G- 1 - G, 1 ) p|0) 

= -2 (0|s 3 (e* p) (G- 1 - G^) (e+p) |o) , s, = 1, s 2 
2ee' 



e 2 + e' 2 



2 ' 



G = H + 1, H = p 2 + VJg), G = p 2 + 1, e' = u-e. (2.2) 



These expressions are similar to the expressions for probabilities of bremsstrahlung 
and pair creation by photon || when the involved charged particles are sub- 
jected to the multiple scattering in a medium. The developed approach is given 
in these papers in detail. The same method was used in for analysis of coher- 
ent scattering of a photon in a medium. In absence of an external field one can 
use the results of this paper putting the momentum transfer A = 0. A modifi- 
cation of radiative correction under simultaneous influence of multiple scattering 
and homogeneous electromagnetic field was considered recently by authors in 
flQfl , where the anomalous magnetic moment of an electron was considered. The 
potential V p (q) can be obtained from Eqs.(2.7)-(2.11) of [ID] with the help of 
standard substitutions to — > — uj,e — ► — e: 



ee' 



V p (g) = -iV(g) + 2k'q, k' = —n 



0J' 



(2.3) 
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where k and F are defined in Eq. (|1.4| ). The potential V(g) has the form 0, |§: 



V{q) = Qg 2 Li + ln— -2C 



Q 



2nZ 2 a 2 es'n a 



m A L) 



Os2 



lSSZ-^e-/, f = f(Za) = {Z a fY. 



\ k(k 2 + (Zaf 



(2.4) 



here C = 0.577216. . . is Euler's constant, n a is defined in ( |1.3|) , g is the coor- 
dinate in the two-dimensional space measured in the Compton wavelength A c , 
which is conjugate to the space of the transverse momentum transfers measured 
in the electron mass m. 



We split the potential V p (g) in the same way as in Eqs.(2.9)-(2.11) of |10 
V p (q) = V pF (q) - iv(g), V pF (g) = -iV c (g) + 2n' g, V c (g) = qg 2 , 



q = QL C } L c = L(g c 



In 



l: s2 



v(g) 



<1Q 



g 



lnf- + 2C . (2.5) 



Ae 2 



\2 n 2> T 
A c6c 

Here the parameter g c is defined by the set of equations (compare with Eq.(2.20) 
of 0): 

q c =1 for 4(k/ 2 + QLx) < 1; 

40c (k' 2 qI + QL C ) = 1 for A{k' 2 + QL X ) > 1, (2.6) 

According to this splitting and taking into account the addition to the potential 

as 



273| ) we present the propagators in Eq.(p^2 



— o 



^pF ~r ^pF ' 



(2.7) 



where 



G p f = T~C p f + 1, G — G p f — iv, H. p f — p 2 + VpF- (2-8) 

The representation of the propagator G~ x permits to carry out its decomposition 
over the "perturbation" v 

G- 1 - G- F = G pF iv G pF + G; F iv G; F ivG~ F + ... (2.9) 

The procedure of matrix elements calculation in this decomposition was formu- 
lated in 0, see Eqs.(2.30), (2.31). Here the basic matrix element is 
< Qi\G~ f \q 2 >■ This matrix element can be obtained from Eq.(2.15) in ||T0 
making the substitutions:^ — > —uj,e — > —e, x/ u k>' 



< Qi\G p f\Qi >= i J dtexp(-it)K pF (g 1 , g 2 ,t), 
o 

K p f(Qi, Q21 1) = K c (g l} g 2l t)K K {g x , g 2 , t), 



(2.10) 



where 

K c (gx,g 2 ,t) 



V 



47ri sinh vt 
K*(Qi, Q 2 ,t) = exp 



iv 



[g\ + g\) coth vt 



sinh vt 



Q1Q2 



Ain /2 t 



v' 



vt 



vt 



2i 



vt 



1 tanh — I K '(f?i + £2) tanh — 



.(2.11) 
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where v = 2y/Tq, q is defined in (|2.5f) . 

In the present paper we restrict ourselves to the main term in the decompo- 
sition (|2.7| ). This means that result will have the logarithmic accuracy over the 
scattering (but not over an external field). The matrix elements entering into 
the mass correction have in the used approximation the following form 



M 1 



1 G p p — Gq 1 1 



1 

47T 



exp(— it) 



v 



sinh vt 



dt. 



M 2 = {0\p(G^-G 1 )p\0 
dt, 



1 °° 
--/exp(-,t) 



r> 



vt 



^(^tanh 2 - 
sinh vt\ v 2 2 



i 



sinh vt y 

M 3 = (0|(e*_p)(G p -i-G 1 )(e + p)|0 



.12 



2nv 



exp(-it) 



<p p = (p(n', v, t) = exp 



Ain'H 





2 , vt 

1 tanh — 

vt 2 



— ^ — tanh 2 — dt, 
sinh vt 2 

(2.12) 



Integrating by parts the term containing 1 / sinh 2 vt together with the sub- 
traction term in expression for M 2 we find 

v °r vt f 4k' 2 \ 

M 2 = / exp(-it)p p tanh — — - + 1 ) dt - M l (2.13) 

Air J 2 \ v 2 J 



Substituting the found expressions for Mi and M 2 into formula ( |2.2| ) and then 
into Eq. (|2.1|) we obtain the general expressions for photon masses squared un- 
der simultaneous influence of multiple scattering in a medium and an external 
electromagnetic field 

2 w J " 

am r oode 



k, 



2tt 



s 2 vip p tajah 



EE' 





vt ( 4k 



r> 



+ 1 - s 3 



sinh vt 



For k\ _ we found respectively 



2am 2 r oodE 



71 



/2 



EE 



1 

7 

vt 



dt 



■tanh 2 — dt 
v sinh vt 2 



(2.14) 



(2.15) 



In the absence of external field (k/ = 0, tp p = 1) we have 



k 2 

K ++ 



2 w J 

am r oudE 



2tt 



EE' 



So^tanh , 

2 2t Vsinh^t t 



dt 



am 2 ( udE 



2tt 



EE' 



1 



Si ( hip - if) ( p + - ) ) + s 2 hp(p) - lnp + — 



(2.16) 
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where p = i/2u, ip{p) is the logarithmic derivative of the gamma function. Sub- 
sisting the result obtained into formula (|1.12 ) we have the probability of pair 
creation which agrees with formula (2.10) in ||. 

In the absence of multiple scattering [y — > 0) we get 



am 



71 



ujde 



EE' 



12. 

S2K t(f F 







dt, 



am 



2tt 



ude ,n 

EE' 



ipptdt, 



(2.17) 



where 



cp F — exp —1- 



(2.18) 



For this case the expressions for k\ 2 (|1.10| ) after substitution results of ( |2.17D 
agree with masses squared of photon calculated in an external electromagnetic 
field (see and references therein). 

We will analyse now the results obtained in different limiting cases. In the 
case when the both characteristic parameters are small iv\ = AQLi < 1, k C 1), 
the main terms of decomposition of the functions k 2 are the sum of independent 
corrections to the photon mass squared both on account of the multiple scattering 
and an external field 



am 



71 



.7 io 



9 uj. 



1 



21Li 



59 
225 



.3\^3n uj 
-1 = — exp 

,2 



8ujf 



3u 



90 



k 



am 



71 



TT UJ 



exp 



Here the notations are introduced: 



m 



+ 



30 



2nZ 2 a 2 X 3 r n a L 1 



U3f 



m- 



-)' 

UJpJ 



rr - m 

e 



uj <ti Ld e , uj uj f .(2.19) 



4.41 • 10 13 Oe. 



(2.20) 



For used notations see Eqs.( |1.4j ), ( p.4| ). The correction ocl/Lx follows from the 
first term of decomposition in Eq. ([2.9|) . In gold the value uj e is uj e = 10.5 TeV, 
this is the typical value for the heavy elements. 

In Fig. the functions Re k+ + (curve 2) and Im k\ + (curve 1) are given for 
the case when the influence of a medium is taken into account only (Eq. fl2.16| )). 
The both curves are normalized to the asymptotics given by Eq. (|2.19|) in the 
limit ujp — ► 00. 

In the case when the influence of an external field is small comparing with 
effect of a medium oj/ojf (1 + uj/uj e ) 3 ^ we can present the expression for k 2 

as 

k 2 = k 2 (uj c ) + Ak 2 (uj F ). (2.21) 
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Here k 2 (uj c ) is the photon mass squared under influence of a medium only (in 
absence of an external field k' — 0) is given by Eq. ( |2.16| ), where the function L c 
is defined in Eqs.( f2T5|) and (|2.6| ) 



Aee' uj 



v = 4iq = i^- — , uj c = uj e — -, L c ~ L x 

UJ Z UJ C L, 



1 / 

1 + — In 1 + - 

2Li \ Co> e 



^2.22) 



Retaining the main terms of decomposition over k' 2 in Eqs. (|2.14j) and ( [2.1 5|) we 
find the corrections Ak 2 (uj F ): 



LU _ OO 



Ak 2 (uj F ) 



am ujuj c 7 de f ( z 

i— 2~ / — / eX P 

27T CJp J UJ J \ V 

^ 



x <s 2 



z / ,z\ 2i 2 z 

tanh - 1 — % — M tanh - 

2 V i// i/ 2 



H ^ — f z — 2 tanh — ) ^ dz, 

z/sinhz \ 2/ 



27T CUp 



z \ tanh 



2 am 2 ujuj c f de 

Ak + _{uj F ) = i— — I —s 3 I exp l . 

uj J \ v) smnz 



dz. 



(2.23) 



o 







In the limit uj u e one can use the asymptotics found in Appendix A in 



for calculation of the photon mass squared k\ 



k 2 ++ {uj c ) 



2am 2 2 .u c 
11(a), a = % — , 

UJ 



7T 



7r 3 a 



11(a) ~ -— + ( 2 In 2 - - ^ In - + A - 

K J 8a V 2) a 48 ' 

1 00 (—l) n 
A = In 2 (2 In 2 - 1) + -(1 + C) - 2 ^ ^ '— Inn = 0.736629. (2.24) 

^ n =l n 



It should be mentioned that in the limit ujf — > oo this formulas g ives at 
uj/uj £ = 10 within accuracy better than 7% and at uj/uj e = 100 within accuracy 
better than 0.3%. 

The corrections Ak 2 (uj F ) due to action of external field in the same limit 
uj j> uj e are 



Ak\Auj F ) 



2am 2 ujuj c 

37T UJ F 



l-ln2), Ak 2 + _(uj F ) 



. am ujuj c 
1 12vr uj 2 f ' 



(2.25) 



In the case when the principal effect is due to an external field, the main 
contribution into the photon mass squared is given by Eq. fl2.17|) and corrections 
connected with influence of multiple scattering are 



k 2 = k 2 (uo F ) + Ak 2 {uj c 



Ak 2 ++ (uj c ) 



am 2 uj r de 



157T UJ r J UJ 

o v o 



exp 



-it I 1 + -k'H 2 



x 



dt + —(2s 2 + s 3 ) \, 



Ak 2 



\UO r 



2am 2 to 

15/T UJ r 



de 

—s 3 

UJ 



exp 







-it I 1 + -K /2 t 2 



2it - 3 



1 



Here the function L c is determined in Eqs. 

2 



1 



3Li 



In 1 + 



and 

UJ \ 



at Q = 



2/t' 2 / 
(2.26) 

(2.27) 



To obtain expressions Eg. ( |2.26[ ) we performed decomposition in formulas fl2.14| ) 
and ( 2.15| ) and integrated by parts the terms with high powers of t in integrals 
over t. Note, that using Eq.( |1.12| ) we can obtain from Eq. ( |2.26| ) the correction to 
the probability of pair photoproduction due to effect of weak multiple scattering. 
The result agrees with Eq.(7.136) in ||. 

In the limiting case u ^> ujp we have from Eq. fl2.17p 



,2 / \ 5am 2 _..T 3 (2/3) f3uj\ 2 / 3 ,. 



[Up 



-A' 



(2.28) 



For corrections originating from the multiple scattering we have from Eq.( [2.26| ) 
the main term of decomposition 



4am 2 _ if T 3 (l/3) /3w\ 1/3 /w F \ 



75tt 



r(2/3) \uj c J \uj c J 



2/3 



Ak 2 + _(u c ) 



(2.29) 



3 Conclusion 



It is curious that in the scope of the used method (see e.g. 0) it is possible to 
find many of obtained in previous section results (up to numerical coefficients) 
basing on very simple form of the amplitude of photon scattering 



M 



a 2 
~ — q 



where q 2 is defined by the equation 



2irujZ 2 a 2 n a L(q 2 + m 2 
q 2 + m 2 



fi 2 

m k 



% 2 



m 



2\2 ' 



(3.1) 



(3.2) 



In the case q 2 -C m 2 {\q\ 2 ~ q 2 ), see Eqs.( |1.3| ) and flO]), the imaginary part 



of the photon scattering amplitude in forward direction (Eqs.( |3~T|) and (|3.2|) )is 
defined by the value q 2 while the real part of this amplitude is a sum of the 
correction ~ Qs/ m2 anc ^ the momentum transfer due to action of a field q 2 F 
(compare with (|2.19| )) 



a 9 
M ~ -m 2 

7T 



m* 



K 

T 



(3.3) 
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At strong multiple scattering (q 2 3> m 2 ) and in the case q 2 F <C q 2 (ft 2 m 6 <C q® ) 
the value q 2 ~ q 2 is defined by Eqs.( |1.3| ) and ( |1.4| ) 

g s %) = 2™Z 2 aV l L(g 2 ), (3.4) 

and the photon scattering amplitude ( |3.1| ) is (compare with Eqs.( |2.24T ) and 



2j| )) 



M ~ — 

7T 



fi 2 

m k 



(3.5) 



When the main effect is the action of the field k ^> m and in the case 
q 2 s {q F ) -C q 2 F (/t 2 m 6 > 9s 9f) the va l ue Qf ~ 9/ we have (see (|L4|)) 

6 2 

gJ(9F) = (3-6) 
and the photon scattering amplitude ( |3.1|) is (compare with Eqs.( |2.28|) and 

(PD) 



M 

7T 



a r 



^q 2 F M + e-^q 2 s (q F ) . (3.7) 



It should be noted that beginning with some photon energy uj = u!^ the 
radiative correction to the value Re Sjk in the absence of a field (n = 0) becomes 
larger than Wq/w 2 (see (|1.1|) ). Let us estimate u>b- According to formulas fl3.1|) 
and ( |3.2| ) this effect originate at values q 2 <C m 2 . Because of this for estimate 
Re M ( |3.1|) it is necessary to take into account the next term in decomposition 
over q 2 in (|3.2j) at k = 0. We find 



Re M s 



Re M s Ti u 2 A 2 



-Im M s 




a; 


7T C<J 


fa 


-^rad 



(3.8) 



For gold one obtains ~ 40 GeV. 

A propagation of high-energy photons in oriented single crystals is one of 
interesting applications of the result obtained above. In this case we have both 
the dense matter with strong effect of multiple scattering and high fields of 
crystal axes or planes. As known, the Landau- Pomeranchuk-Migdal (LPM) effect 
(influence of multiple scattering on processes of bremsstrahlung and pair creation 
by a photon) is most pronounced in the heavy elements. The same is valid for 
the process under consideration. Let the high-energy photon incident on crystal. 
The angle of incidence is small and such that the distance from axis g (or the 
distance from plane x) can be considered as a constant on the formation length 
of process (the constant field approximation is applicable, see Sections 12,15 in 

I)- 

For orientation of a crystal along an axis the ratio of density of atoms in the 
vicinity of axis n(g) to the mean density n a is 

n a 2nu 2 dn a ' 
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where U\ is the amplitude of thermal vibrations of atoms, d is the mean distance 
between atoms which form the axis. This ratio is maximal at q = 0. For 
numerical estimates we use for definiteness the tungsten single crystal. For the 
axis < 111 > in W the ratio £ aa; (0)=370 at the room temperature (T = 293 K) 
and £ aa; (0)=1020 at T = 77 K. The effect of multiple scattering becomes strong 
at characteristic photon energy u e (n a ) ~ 11 TeV and this value is inversely 
proportional to the density. So we have that uj e (g = 0) ~ 30 GeV at T = 293 K 
and uJ e {Q = 0) ~ 11 GeV at T = 77 K. It should be noted that within logarithmic 
accuracy we neglect by relatively small variation of L\ due to substitution the 
screening radius a 2 s by the value 2u\. 

It is useful to compare these estimates with known threshold energies tut at 
which the probability of pair creation in the field of axis is equal to the probability 
of the Bethe- Heftier mechanism, see Table 12.1 in ||. For photon energy uj > uj t 
the process of pair creation in the field of axis dominates. In W crystal for 
< 111 > axis oj t = 22 GeV at T=293 K and oj t = 13 GeV at T=77 K. It 
is seen that for these energies the ratio u)/u) e which characterize the strength 
the LPM effect is of the order of unity. At uj ~ uj t the maximal value of the 
parameter k(q) which determines the probability of pair creation by a photon in 
a field is also of the order of unity (at g ~ U\, see the mentioned Table). So we 
reach the conclusion that at some energy (for axial orientation of crystal) all the 
discussed effects are essential simultaneously. The analysis in this situation will 
be published elsewhere. For example, to calculate the influence of the field of axis 
on the polarization tensor one have to average the general formula ( p. 14 ) over 



all values of g (this is integration over d 2 g with the weight n±, where n± = n a d 
is the density of axis in the perpendicular to them plane). 

At planar orientation of crystal the ratio of the density of atoms in the vicinity 
of plane n(x) to the mean density n a is 

n a v27TMi 



where d v \ is the distance between planes. For the plane (110) in W crystal at 
T = 293 K one has £ p z(0) ~ 18 and effective value uj e ~ 600 GeV. This value 
is ~ 2.5 higher than the threshold energy u t for this plane (see Table 15.1 in 
0). Let us note that in both axial and planar cases we made estimations for the 
maximal value of £(0). 

In the crystals where the atomic number Z is not very high (Ge, Si, C) 
the ratio co> t /c<j e (0) is smaller than unity. So for uj > uj t one can use Eqs. fl2.17j ) 
and (|2.26|) while for u uj t the formula (|2.19|) is applicable in which along with 



known results (see e.g. 0) there is the new term in Re k'] r+ which is proportional 

to (u/u e ) 2 . 
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Figure caption 



The functions Re k\ + (curve 2) and Im k\ + (curve 1) versus the photon energy 
taken in units u e (because of this the curves are universal) for the case when the 
influence of a medium is taken into account only (Eq.( [2.16| )). The both curves 
are normalized to the asymptotics given by Eq. fl2.19| ) in the limit up — > oo. 
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